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, The contraction processes of metal-free molecular clouds of starlike mass (or cloud cores) 

• are investigated. We calculate radiative transfer of the H2 lines and examine quasi-static 

contraction with radiative cooling. Comparing two time-scales, the free-fall time iff and the 
time-scale of quasi-static contraction t qsc (~ t C ool, the cooling time) of these cores, we find 
that the ratio of the two time-scales tff/tqsc, i.e., the efficiency of cooling, becomes larger 
with contraction even under the existence of cold and opaque envelopes. In particular, for 
fragments of primordial filamentary clouds, for which ts ~ t qsc at the fragmentation epoch, 
they collapse dynamically in the free-fall time-scale. This efficiency of cooling is unique to 
. line cooling. 

m . 

(N ■ 

§1. Introduction 

00 : 

One of the most conspicuous features of galaxies is the existence of stars. To 
discuss the formation and evolution of galaxies, then, it is indispensable to investi- 
Qh' gate star formation processes in protogalactic clouds. In addition, stars play crucial 

roles in galactic activities, for example, ultraviolet radiation and supernova explo- 
sions, which are important in some cosmological contexts. Recently, considering the 
hierarchical clustering scenario, several authors stressed the importance of stars in 
pregalctic objects for the reionization of the universe (e.g., Ref. ||)). However, the 
.,— 1 \ formation of first stars is one of the most poorly understood processes involved in 

galaxy formation. For this reason, theoretical predictions on the reionization by first 
stars suffer great uncertainty. Therefore it is definitely well-timed to investigate the 
formation process of first stars. 

The evolution of primordial clouds from which first stars form has been investi- 
gated with a one-zone approximation by many authors, a However, the mass scale 
of the clouds which are able to collapse is much larger than the typical stellar mass. 
Then, we must discuss the fragmentation processes of protogalactic clouds and the 
evolution of these fragments into stars to understand the star formation. 

With respect to the fragmentation process, several authors have envisaged that 
protogalactic clouds first collapse disk- like, and then fragment into filaments (e.g., 
Ref. |3|) ) . Uehara et al. tf investigated the collapse of filamentary primordial clouds 
with a one-zone approximation and concluded that the minimum mass of the frag- 
ments becomes ~ 1M©, that is, essentially the Chandrasekhar mass. Their calcu- 
lation also showed that, for low mass fragments, almost all of hydrogen atoms are 
converted into molecular form before the fragmentation occurs. 
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On the other hand, Stahler et al-H* studied the main accretion phase of a pri- 
mordial protostar assuming a very low mass hydrostatic core (i.e., a stellar core, see 
Ref. |6|)) initially, and stationary accretion of the envelope onto it. They adopted a 
higher mass accretion rate (4.41 x 10 _3 Mq yr _1 ) than the present-day counterpart 
(1O" 5 M yr" 1 , e.g., Ref. 0)). This higher mass accretion rate is caused by higher 
temperature ~ 10 3 K of primordial gas clouds. However, the question of how a frag- 
ment of star-like mass reaches this phase is not considered and still remains unclear. 
In this paper, we study the evolution of stellar mass primordial clouds shortly after 
the fragmentation. 

How primordial clouds of stellar mass, which are fragments of filamentary clouds, 
evolve after the fragmentation depends on their efficiency of cooling, i.e., the ratio of 
the free-fall time-scale tg to the cooling time-scale t coo \ (e.g., Ref. If % < i coo i 
(i.e., cooling is not effective), while the cloud cools, it has enough time to adjust 
itself to a new hydrostatic configuration. Then it contracts, maintaining a nearly 
hydrostatic equilibrium (i.e., Kelvin-Helmholtz contraction). On the other hand, in 
the case ifj > t C ooh the cloud collapses dynamically in its free-fall time-scale after 
substantial cooling, and no hydrostatic equilibrium is established. In the latter case, 
the cloud may fragment into smaller pieces again. 

At the time of the fragmentation of filamentary clouds, the two time-scales of 
fragments, t coo i and tg, are of the same order of magnitude,^ since a filamentary 
cloud fragments when the time-scale of contraction idyn = plidp/dt) becomes nearly 
equal to the time-scale of fragmentation if rag = 2A/^/2^rGp,w where p is the density 
of the cloud. At the epoch of fragmentation, the contraction proceeds isothermally 
with regard to time. Then td yn ~ (7 a d — l)t C ooi> where the adiabatic exponent 7 a d 
is 5/3 for monatomic molecules and 7/5 for diatomic molecules. Since the fragmen- 
tation time-scale if rag of a filamentary cloud is identical (up to numerical factors) 
to the free-fall time-scale of a spherical fragment, t coo i ~ tg for the fragment at 
the fragmentation epoch. Therefore, the nature of the evolution depends on whether 
the ratio £ C ooiAff [more precisely, t qsc /ts, where t qsc is the time-scale of quasi-static 
contraction (see §3 below)] becomes larger or smaller as a fragment contracts. 

To estimate the cooling time, we must investigate the cooling processes in the 
primordial clouds. The principal difference of a primordial star forming region from 
its present-day counterpart is its lack of heavy elements and consequently dust grains, 
which are efficient coolants, keeping present-day molecular clouds almost isothermal 
(~ 10 K) over many orders of magnitude in density (e.g., Ref. [l0|)). Therefore, when 
we investigate primordial star formation, it is no longer possible to impose isother- 
mality. We must study cooling/heating processes and solve the energy equation to 
determine the temperature. At low temperature (T < 10 4 K), the Lya line of the 
neutral hydrogen is hardly excited, and there is no atomic cooling process. If there 
are H2 molecules, however, these clouds can cool by H2 rotational/vibrational tran- 
sition lines. According to Uehara et al.,0' at the time of fragmentation, filamentary 
clouds are optically thick for line center frequency of some lines. This should also 
hold for fragments of them. The surface temperature of the clouds is low, in com- 
parison with the central or averaged temperature. Consequently, the efficiency of 
cooling by optically thick lines may decrease. For example, the continuum energy 
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flux decreases as T% 



where T e g is the effective (photospheric) temperature of a 



cloud. In an optically thick case, the photosphere is located near the surface, and 
the low surface temperature causes a severe decrease in the cooling rate. However, in 
the case that cooling is dominated by optically thick lines, optical depth is not only 
different from one line to another line, but it also varies with frequency even within 
one line. In fact, some lines are optically thick at their line center frequencies, but 
clouds are totally transparent to intervening frequencies between lines. Therefore 
the notion of a photosphere cannot be easily implimented. Final processed line pro- 
files depend on the structure of the cloud and cannot be known before we calculate 
radiative transfer resolving frequencies within a single line. Moreover, the number 
of effective lines changes as the cloud contracts. For these reasons, to understand 
the evolution of the fragments, it is necessary to treat the H2 line transfer problem 
in detail. 

The outline of this paper is as follows. In §2, we describe our scheme of the cal- 
culation of the luminosity and the cooling rate via the H2 emission lines for spherical 
clouds and give some results. In §3, we compare the time-scales ts and t coo i for gen- 
eral initial conditions and discuss the evolution of fragments of filamentary clouds. 
Some interpretations of the results are also presented in §4. Finally, we summarize 
our work and discuss its implications with respect to the primordial star formation 



§2. Radiative processes in primordial molecular cloud cores 

We consider spherically symmetric gas clouds of starlike mass (~ 1M@) ('molecu- 
lar cores') composed only of H2 molecules, because the hydrogen atoms are converted 
into molecular form before the fragmentation for these low mass fragments, o 1 We 
can neglect helium, since it is thermally inert at these low temperatures. As men- 
tioned above, these clouds lose their thermal energy by radiative cooling via the H2 
rotovibrational lines, and presumably contract into protostars. 

In this section, we describe our calculation scheme for the luminosity and cooling 
rate of a spherically symmetric cloud. The specific intensity I v (ergs sec _1 cm _2 sr _1 Hz _1 ) 
along a ray is calculated by solving the radiative transfer equation (e.g., Ref. pi])), 



Here s is the displacement along the ray, and a v and j v are the absorption and 
emission coefficients, respectively. These coefficients can be written by using Einstein 
A— and B— coefficients: 



in §5. 



dh 
ds 



V 



(2-1) 






Here m and 112 are the number densities of the molecules in the lower and upper 
level of the transition, respectively, and <j>(y) is the line profile function. In our 
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calculation, the line center optical depth of the lines across a cloud r c is no more 
than several hundreds. Since Lorentz wings of a line become more important than 
its Doppler core only when r c is as large as 10 3 , we only consider line broadening 
owing to the thermal Doppler effect, i.e., 

Mv) = - 1 e -("-"o) 2 /(^D) 2 . (2.4) 
Here vq is the line center frequency, and the Doppler width Av-q is defined by 



vo 2kT . . 

Avd = —\ , 2-5 

c y /iron 

where T is the temperature at each radius in the cloud, tbh = 1-67 x 10~ 24 g is 
the mass of a hydrogen atom, and fi = 2 is the molecular weight of the molecular 
hydrogen. We include the first three vibrational states, with rotational levels up to 
j = 20 in each state, following Palla et al.cP We take ^4-coefficients of molecular 
hydrogen from Turner et al. 0) -B-coefficients can be obtained from A-coefficients by 
the following Einstein relations, 

31-612 = 32-821, (2-6) 

A 2 i = —§-B2u 2-7 
<r 

where gt is the statistical weight of the level i. The critical number density n cr , where 
radiative and collisional deexicitation rates become equal for molecular hydrogen, is 
~ 10 4 cm -3 . For a typical number density of fragments (~ 10 10 cm" 3 3> n cr ), almost 
all excited molecules would be deexcited by collision with other molecules. Therefore, 
collisional local thermodynamic equilibrium (LTE) is established, and scattering of 
photons can be neglected. These properties make the transfer problem of H2 lines 
fairly tractable. 

Figure 1 displays the geometry of grids of the radiative transfer calculation in 
the case of spherical symmetry (Ref . |l3|) ; see also Ref . |l4|)). To know the luminosity 
at the z-th radial grid rj, we solve the transfer equation ( |2-1[ ) with frequency v along 
from the 0-th to (i — l)-th rays, and we obtain the specific intensity at rj in the 
direction 0^, I„(ri,9ij) {—i + 1 < j < i — 1). This procedure will be repeated for 
each frequency mesh. We consider the frequency range for each line within four times 
the Doppler width corresponding to the central temperature from the line center. 
Each line is usually separated into 20 frequency meshes (in some cases, e.g., Fig. 7, 80 
frequency meshes), and the number of radial meshes n = 50 in our runs. Integrating 
over and v , we obtain the monochromatic energy flux F v and the luminosity L: 



F v (r) = 2tt [" I v (r,O)coB9d0, (2-8) 
Jo 

roc 

L(r) = 4vrr 2 / F v {r)dv. (2-9) 
Jo 
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As for cooling rate per unit mass A(m), we differentiate the luminosity L with 
respect to the Lagrangian mass coordinate m, 

OL , 
A(m) = —, (2-10) 



m (r) = [ r Aitr' 2 pdr'. (2-11) 
Jo 



Figure 2 displays the contribution from each line to the luminosity at the surface 
of a typical cloud. The mass of the cloud is 1M , the number density at half mass 
radius is 10 11 cm -3 , and the density /temperature distribution is represented by the 
Emden function of polytropic index 2.5. In this spectrum, there is a gap around the 
wavelength of 4/xm, and lines can be classified into longer and shorter wavelength 
groups. The shorter and longer wavelength lines correspond to transitions with and 
without changes of vibrational levels. 



§3. Evolution of fragments 



As initial conditions, we assume polytropic gas spheres in hydrostatic equilib- 
rium; namely the density/temperature distribution is represented by the Emden 
function of polytropic index N. Clouds are cut off at the radius r s , where the 
density falls off by a factor of 10 -3 from the central value. Therefore, parameters 
characterizing initial conditions are the effective polytropic index N, the total mass 
of the cloud M, and the number density at the half mass radius. For a given 
initial configuration, we can obtain the specific entropy distribution s(m,t = 0). 

We then bring forward the cloud one time step At quasi- statically in the following 
way. We calculate the cooling rate A(m) as described in §2, and advance the specific 
entropy distribution s(m, t) by one time step At using the heat equation 

With the entropy distribution at the next time step s(m, t + At) = s(m, t) — 
(A(m)/T)At, we find the new equilibrium configuration using the equations of hy- 
drostatic equilibrium, 

dr 1 
dm 4nr 2 p ' 

dp Gm 
dm 4irr 4 ' 

and the equation of state, 

p cc exp[s/c v ]p 7ad . (3-4) 

Here 7 a( j = 7/5 is the adiabatic coefficient and c v = - d 1 _ 1 — ^ is specific heat of hy- 
drogen molecule at constant volume. We impose a constant external pressure at the 
outer boundary. For dense and cold clouds, like the fragments under consideration, 
radiation pressure can be neglected, and the pressure can be assumed to be given by 
the gas pressure only. 



(3-2) 
(3-3) 
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For given states of clouds, we compare two time-scales, i.e., the free-fall time 



/ 3tt 

and the time-scale of quasi-static contraction, 

VM = pI (%) ■ (3-6) 

\OZ J m, quasi-static 

Here, p{m) is the mean density within the mass coordinate m, and each quantity is 
evaluated at the coordinate m. From the virial theorem, t qsc is of the same order of 
magnitude of t coo \. 

Figure 3 shows the contours of the ratio tg / t qsc at the center (m = 0) of clouds 
with (a) N = 2.5 (adiabatic stratification) and (b) N = 8. Regions around the 
curves tg/t qsc = 1 are expected to be initial states of fragments. Although Fig. 3 
displays the values at the center, this ratio does not change significantly within a 
cloud, except near the surface. The upper-right regions correspond to higher central 
temperature > 2000 K, where the molecular hydrogen begins to dissociate. From 
the two panels of Fig. 3, we can see that the ratio ts/t qsc depends only slightly 
on the effective polytropic indices N, rather than on masses and densities. Note 
that for clouds of the same mass, the ratio tg/t qsc is greater for denser cloud. This 
implies that as a cloud contracts, the ratio tg/t qsc becomes larger. In particular, 
for fragments of filamentary gas clouds, for which tg ~ tq SC initially, tqsc becomes 
shorter than tg as they contract; i.e., such clouds collapse dynamically. Note that 
when t qsc < iff, clouds do not contract quasi-statically, and i qsc does not possess the 
meaning of the collapse time-scale. In this case, t qsc merely measures the time-scale 
of cooling. 

The discussion above is for fixed N. However, these clouds change their entropy 
distribution as they contract. Therefore it is not enough to discuss their evolution 
only with fixed N. We need to calculate the evolution of clouds contracting quasi- 
statically and to show that they in fact begin to collapse dynamically. We choose 
initial conditions under which clouds contract quasi-statically and pursue their evo- 
lution. Their initial parameters are (a) M = O.5M , = 10 9 cm -3 , N = 2.5 (the 
asterisk * in Fig.3(a)), and (b) M = 0.2M Q , n h = 10 9 cm™ 3 , N = 20. For the 
parameter representing the structure of clouds, we employ the effective exponents 
for structure, 

_ (dlnp/dm) t 

7str(m) = — 7^-. 3-7 

(dmp/dm)t 

Figure 4 displays changes of the 7str("i') as the clouds contract. We can see 
that values of 7 s tr("7') become nearly flat and approach 1.12 — 1.13 when the central 
temperature reaches ~ 2000 K, which is the dissociation temperature of H2. These 
values of 7 s t r correspond to N ~ 8, regardless of the initial N values. However, 
for the fragments, the clouds start to collapse dynamically before these states are 
reached. In such a case, this limiting value of N has no significant meaning. But 
the very existence of the limiting value supports our discussion above with fixed N. 
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In fact, from Fig. 5, where the change of time-scales and their ratio are shown for 
the case (a) of Fig. 4, we can see that the cloud indeed collapses dynamically after 
it contracts quasi-statically to some extent. 

Although 7 s t r becomes ~ 1.13, the effective exponent for temporal evolution of 
a fixed mass element, 

(dlnp/dt) m 



7cvO) = 



(3-8) 



(<91np/<9f) m ' 

has different value. The distributions of 7 CV at different times are shown in Fig. 
6. Effective adiabatic coefficients, 7 CV , are always nearly equal to 4/3 in the case of 
quasi-static contraction. The reason for this is that the relation between 7 CV and 
7ad (= 7/5) is 



7e 



7ad 



'■qsc 
^cool 



(3-9) 



where i coo i = e/(dL/dm) and e(= - d-1 is internal energy per unit mass. From 
the virial theorem, 

3(7ad - + E e = 0, (3-10) 

we find that 



F/ dEi 



dE e 



dt 



where E\ and E g are the total internal energy and gravitational energy: 



rM r 

E{= edm, E g = — 
Jo Jo 



M Gm 



dm. 



(3-H) 



(3-12) 



In the case of homologous contraction, Eg oc p 1 / 3 , where we compare p for a fixed 
mass element (a homologous point), and t qsc becomes constant along the radial 
coordinate: 



t. 



qsc 



On the other hand, by using total energy of a cloud, 

'4 



W = E { + E g 



7ad E u 



its total luminosity can be given as 

T dW 
dt 

The mean cooling time-scale becomes 

(t ) Ei 

\ L coo\/ 





dEi 


(7ad " |) 


dt 



E/ 



dEi 
~~dt 



L 3( 7ad - |) 

Therefore, for homologously contracting spherical clouds, we have 

tqsc — I 7ad _ g ) V'cooU • 



(3-14) 



(3-15) 



(3T6) 



(3-17) 
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For the case of non-homogeneous contraction, this equation is not exactly valid. 



However, holds approximately in this case. From Eq.( |3-9|) , ( 3-17 ) and identifying 
global cooling time-scale (t coo \) with local one t coo \, j ev becomes nearly equal to | 
as a whole. From Fig. 6, it can be seen that the initial steep gradients of the 7 ev are 
levelled, and they assume values around 4/3. This is because there exists a limiting 
structure for clouds contracting quasi-statically, and then the contraction becomes 
roughly homologous. 



§4. Interpretation of results 



In §3 above, we have shown that the stellar mass molecular cores collapse dy- 
namically, although some lines are optically thick at line center. This fact indicates 
that the efficiency of cooling by H2 lines does not fall significantly even in such a 
situation, in contrast to the case of continuous radiation. The reason for this can be 
interpreted as follows. 

Figure 7 displays the change of a processed profile of a single optically thick 
line (optical depth at line center across the cloud is 76) along a ray passing though 
the center of the cloud. Since the cloud is optically thick at line center, the specific 
intensity at line center is saturated with the blackbody value determined from the 
local temperature. On the other hand, optical depth at wing is smaller than unity. 
The final processed profile at the cloud surface is double peaked, and the peak 
frequencies correspond to frequencies where the optical depth across the core r VjCor e > 
1 (i.e., the radiation is saturated with the value of the core temperature) and that of 
envelope T Ufinv < 1 (i.e., the radiation can be larger than the value of the envelope 
temperature). The height (i.e., the normalization of intensity) and width of these 
peaks are of the same order of magnitude as the blackbody value determined from 
the temperature near the center, not around the surface of the cloud. The reduction 
factors from the central values are only several. Owing to the effect described above, 
which is a distinctive feature of line cooling, the temperature around the cloud 
center can be 'seen' from outside the cloud, even when lines are optically thick at 
line center and the surface temperature is relatively low. Accordingly, its cooling 
proceeds efficiently, as long as there are enough frequency ranges where r VjCor o > 1 
and r^env < 1- Therefore, the decrease in the efficiency of cooling caused by the low 
surface temperature is not severe in the line cooling case, in contrast to the case of 
cooling by continuum. 

Figure 8 displays the ratio of Li, the contribution to the luminosity at the surface 
from each line i, to Lthin.i; the value in the "optically thin" case (i.e., L{r s ) = 
Ei= a ii lines L h and L mn>i = 4vr/ wholc cloud dV J hnc ^vjy), versus optical depths at 
line center across the clouds r Citot for each line. All marks of a given variety represent 
lines of a single given cloud, and the solid curve represents the value of Lj/£thin,i 
for the uniform cloud of 1M , whose number density equals 10 11 cm -3 . From this 
figure, it can be seen that Lj/L t hi n ,i falls only as ~ T tot _1 . This simply reflects the 
fact that in the "optically thin" case, the specific intensity along a ray at optical 
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depth r (even if r > 1) is determined as 

/,,thin= f T S v (T')dr', (4-1) 







where S v = j u /& u is the source function. On the other hand, in the optically thick 

case, 

h ~ S u (t). (4-2) 

Then, if r > 1 

Within ~ t~\ (4-3) 

Therefore, no significant suppression of cooling caused by the low surface tempera- 
ture arises in Fig. 8. 

As discussed above, the effects of the surface temperature are not severe even 
for optically thick lines. Thus it is appropriate to use the value of a homogeneous 
cloud in estimating the cooling time. This is given by 

Ei _ MCyT M l/3 T -3 2/3 f - 1 / 4 . 41 

where Ei, M,T, R,c v and n are the internal energy, the mass, the temperature, the 
radius, the specific heat at constant volume, and the number density of the homo- 
geneous cloud, respectively. The quantity / is the ratio of the luminosity of the 
cloud to the blackbody luminosity of the temperature T. In the case that cooling 
is dominated by line cooling, by using the line width Av and the effective number 
of lines a c , f can be written as / = L j '4tt R 2 aT A ~ a c 4r- We adopt the Doppler 

width = ^-y ^lroH ^ or ^ e ^ ne w i°-th- Assuming the cloud is almost in hydrostatic 
equilibrium initially, we can relate the temperature with the mass and the number 
density: T oc M 2 ^n 1 / 3 . From these relations, we obtain 

i coo i oc M~ 2 ri~ 1/2 a~ l . (4-5) 

On the other hand, the free-fall time is given by 



Then the ratio of these two time-scales is 

oc M' A a c . (4-7) 



^cool 

As the cloud contracts, the number of optically thick lines (in this case, a c ) increases, 
and the ratio ts/t coo \ becomes larger. This accounts for the results of numerical 
calculation described above (see also Ref. and 15)). 
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§5. Conclusion and Discussion 

Our principal result is that if cooling is dominated by H2 line cooling, a primor- 
dial cloud of stellar mass scale is able to cool efficiently. This is because the efficiency 
of line cooling does not decrease significantly by the effect of a line profile, even when 
the cloud is optically thick for line center frequency and its surface temperature is 
relatively low. In particular, fragments of filamentary clouds collapse dynamically 
after fragmentation. Although our calculation is limited to ~ l-M© composed only 
of H2, our results can be applied to more general molecular cores, as long as cooling 
is dominated by optically thick line emissions. 

A dynamically collapsing spherical gas cloud is unstable to a non-spherical per- 
turbation and becomes disk- like. 113 Therefore, dynamically collapsing fragments may 
also become disk-like. Since a self-gravitating disk is unstable to fragmentation, 
cloud cores may fragment again. Whether repeated fragmentations are possible or 
not depends on the thickness of the disk. Fragmentation occurs for a thin disk (axis 
ratio > 2n), but does not occur for a thick disk.HZP Because fragments of filamentary 
clouds start collapsing from nearly hydrostatic equilibrium, it seems unlikely that 
they reach a highly flattened configuration. Therefore it is expected that even if 
fragmentation were to occur, it would not be the case that numerous small pieces 
would be produced. 

Uehara et al.& investigated the fragmentation of filamentary clouds with a one- 
zone approximation. Recently, Nakamura and UmemurallM' examined cylindrical 
collapse of primordial clouds with radial structures and showed that at the epoch 
of fragmentation the clouds display clear core-envelope structures. They suggested 
that Uehara et al.'scP one-zone clouds correspond to the cores of their clouds. Hence, 
if the contraction of fragments were slow, the cores would accrete ambient gas. This 
would result in more massive ones. However, our analysis shows the contraction 
is sufficiently rapid (its time-scale is the free-fall time-scale in the core). Thus the 
accretion of the envelope (its time-scale is the free-fall time-scale in the envelope), 
if it exists, becomes significant only after the molecular core contracts sufficiently 
into a stellar core. The accretion of the envelope does not affect the dynamics of the 
contracting core significantly. Consequently, it is possible to treat the contraction of 
the core and the accretion of ambient gas separately. 

In this work, we showed that fragments of filamentary clouds collapse dynam- 
ically, although we have not pursued yet the actual evolution of the cloud after 
they begin dynamical collapse. To elucidate the actual contraction of primordial 
protostellar clouds, elaborate hydro dynamical calculations are needed. This will be 
presented in a forthcoming paper. 113* 
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Fig. 1. Ray geometry of radiative transfer calculation in spherical symmetry. We solve frequency- 
by-frequency radiative transfer equations along rays tangent to each mass shell and obtain the 
intensity I v {ff) at each mesh point. 



Fig. 2. The contribution from each line to the luminosity at the surface versus its wavelength. 
The mass of the cloud is IMq, the number density at half mass radius is 10 11 cm~ 3 , and the 
density /temperature distribution is represented by the Emden function of polytropic index 2.5. 



Fig. 3. Contours of the ratio ts/t qsc at the center of clouds of polytropic index (a) N = 2.5 and (b) 
N = 8. The contour spacing is logarithmic with an increment of 0.5. The upper-right regions 
correspond to the temperature T > 2000 K, where hydrogen molecules begin to dissociate. 



Fig. 4. Evolution of the effective exponent for structure, j B ti(m) = (g^p/g^)' ' ^ or ^ wo c l uas ^ 
statically contracting clouds. The initial conditions (solid lines) are (a) M — O.5M , n h — 
10 9 cm -3 , N = 2.5. (b) M = 0.2M©, n h = 10 9 cm" 3 , N = 20. Each curve depicts 7 str (m) at 
some specific time (in order of time: solid, dotted, short-dashed, long-dashed and dash-dotted 
line). The central temperature and density are (a) initial: 149 K, 1.35 x 10~ 14 gem -3 , final: 
2048 K, 1.50 x 10" 10 gcm~ 3 and (b) initial: 89 K,2.48 x 10~ 14 gem -3 , final: 2154 K.4.41 x 
10~ 9 gem -3 . These values approach 1.12 — 1.13 around the central region for both clouds when 
the central temperature reaches ~ 2000 K. Note that both clouds start to dynamical collapse 
before reaching the final states shown in the figure. Therefore, our quasi-static formulation 
cannot be applied, and thereafter the evolution shown in the figure does not correspond to the 
actual evolution of the clouds. 



Fig. 5. Evolution of free-fall time ts and dynamical time for quasi-static contraction t qBC at the 
center and their ratio for a cloud which contracts quasi-statically at the beginning (case (a) of 
Fig. 4). The time-scales tg and t qBC and the time ellapsed from the begining t are all normalized by 
iqsc.o, the initial value of t qac - This figure indeed shows the cloud begins to collapse dynamically 
after ~ 1.5i qsc ,o- 



Fig. 6. Evolution of the effective exponent for temporal evolution, "y cv (m) = ^q^J/q^" 1 , for the 
same clouds as Fig. 4. Here the same line corresponds to the same time as Fig.4. 



Fig. 7. Processed line profiles of an optically thick line along a ray passing through the center of the 
cloud at the radii containing % (center, solid line), 1.4 % (dotted line), 9.9 % (short dashed 
line), and 100 % (surface, long dashed line) of the total mass. The optical depth at the line- 
center frequency r c = 38, 46, 54and76, respectively, measured from the surface of incidence. The 
line corresponds to the transition from (u = 1, j = 3) to (i> = 0, j = 1). The parameters of the 
cloud are the same as in Fig. 2. The abscissa represents frequency from line center normalized 
by the Doppler width at the center of the cloud Avo,c- 
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Fig. 8. The ratio of Li, the contribution to the luminosity at the surface from each line i, to Lthin.i, 
the value in the "optically thin" case, as a function of the optical depth at line-center frequency 
across the clouds r c ,tot for each line. These relations are demonstrated for two clouds. Open 
circles and crosses correspond to lines from the cloud of (N = 2.5, M — 1M©, nu — 10 11 cm" 3 ) 
and (N = 8, M = IMq, n h — 10 11 cm" 3 ), respectively. The solid curve represents Li/I/thin, i 
for uniform clouds. 
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